Ultra-high energy cosmic rays, cosmological constant, and 6- vacua 
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We propose that the origin of ultra-high energy cosmic rays beyond the GZK cutoff and the 
origin of small cosmological constant can both be explained by vacuum tunneling effects in a theory 
with degenerate vacua and fermionic doublets. By considering the possibility of tunneling from 
a particular winding number state, accompanied by violation of some global quantum number of 
fermions, the small value of the vacuum dark energy and the production of high energy cosmic 
rays are shown to be related. We predict that the energy of such cosmic rays should be at least 
5 x 10 14 GeV. 



The two outstanding puzzles of modern astro-particle 
physics are the observed value of the small cosmological 
constant, and the origin of ultra-high energy cosmic rays. 
The latter with energies ranging from 10 5 — 10 22 eV, ever 
since it was observed in the first half of the last century, 
posed an open question which attracted many new ideas 
within conventional astrophysics; from the particle spec- 
trum of the Standard Model (SM) to beyond the SM [1]. 
Undoubtedly such a vast range of energies could never 
be covered by a single source for the origin of the cosmic 
rays. The observed broken power law spectrum of cos- 
mic rays gradually steepens as the energy increases from 
10 15 eV, known as a knee, to 10 18 eV known as an ankle, 
and subsequently flattens above 10 18 eV. It is usually be- 
lieved that the first transition in the observed spectrum 
reassigns the origin of cosmic rays from galactic to extra 
galactic sources. However there are obvious constraints 
on a primary particle accelerated up to 10 18 eV if they 
are either charged or a heavy nuclei which interacts with 
a cosmic micro wave photon background T ~ 3 K, and 
thus they cannot traverse further than a few Mpc with- 
out losing energy. This is known as Greisen-Zatsepin- 
Kuzmin (GZK) effect [2,3]. Therefore distant astrophys- 
ical sources which might be able to generate such ener- 
getic particles might not be a suitable candidate for ultra 
high energy cosmic rays with energies beyond 10 18 eV. A 
simple solution to this impasse is to look for a candidate 
which is not only capable of producing ultra high energy 
cosmic rays but also avoids the GZK cut-off. 

Such candidates could be topological defects [4], or 
they could come from the decays of the primordial super 
heavy particles [5] . In the latter scenario the mass of the 
unknown X particles could be ranging from 10 12 GeV 
and above. However the longevity of the X particles 
(equivalent to the age of the Universe ~ tx ~ 10 10 years) 
demands an extraordinary suppression in their interac- 
tion. Rather interesting solutions have been put for- 
ward [6,7]. In the latter reference it was assumed that the 
required life time can be obtained by imposing discrete 
gauge symmetries even if X is an elementary particle. 
Similar ideas in string theory can be found in [8] . While 
in [6] , the reason for the long life time of the X particles 
was explained via the instanton mediated decay, which 
we shall explore here in some detail. 



An interesting connection can be made between the 
abundance of the non-luminous cold dark matter, ob- 
served as 30% of the critical energy density of the Uni- 
verse: p c 4 x 10 -47 (GeV) 4 , and the origin of the 
cosmic rays (above GZK cut-off) provided that the cold 
dark matter constituent is X particles with a mass mx ~ 
10 11 — 10 15 GeV. Such heavy particles can be produced 
abundantly to match the correct cold dark matter abun- 
dance right after the end of inflation [9], or from the 
direct decay of the inflaton [10]. 

On the other hand the majority of the energy den- 
sity ~ 70% is in the form of dark energy, whose con- 
stituent is largely unknown, but usually believed to be 
the cosmological constant [11]. The bare and the ob- 
served cosmological constant is a severe problem, espe- 
cially why the observed cosmological constant is so small 
~ 2.8 x 10~ 47 (GeV) 4 , or in other words in Planck units 
~ 10~ 120 (M p ) 4 (where we use the reduced Planck mass 
M P ~ 2.4 x 10 18 GeV). One would naively expect that 
even if the bare cosmological constant can be made to 
be vanishing, the quantum loop corrections would even- 
tually lead to quadratic divergences ~ M 4 . In other 
words why and what keeps the cosmological constant so 
small as we see today? Despite many attempts in a con- 
ventional big bang cosmological set up [12], a convinc- 
ing solution is still elusive. Recently it was pointed out 
in [13] that, the present observable cosmological constant 
can be obtained if the original vacuum can be associated 
with a non-trivial winding number. It was also assumed 
that the bare cosmological constant and the vacuum en- 
ergy density vanished not in any one specific vacuum 
but in the superposed, or 0-vacuum at some high scale 
- 4 x 10 16 GeV. 

In this paper we address some interesting cosmological 
consequences of the 0-vacuum. We propose that if the 
X particles decay can be explained via instanton media- 
tion due to the transition from one non-trivial vacuum to 
another, then the superposed vacua can also be responsi- 
ble for generating a non- vanishing but small cosmological 
constant as we observe now. As a result, we can relate 
the origin of the cosmological constant with the ultra- 
high energy cosmic rays. Regarding the 9- vacuum we 
also generalize the description of Ref. [13], where it was 
strictly assumed that the state of the Universe is solely 
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given by the \n = 0) vacuum. In what follows, we as- 
sume that the quantum state of the Universe can be a 
superposition of the n- vacua of the non-trivial SU(N). 
As we shall see we will have a concrete prediction for the 
mass of the fermion, which may act as a source for the 
ultra-high energy cosmic rays. We begin by considering 
the origin of the 9- vacuum. 

Apart from the familiar invariance principle of gauge 
theories under small gauge transformations (those con- 
nected continuously to the identity) , it is well-known that 
most grand unified theories such as SU(N) or SO(N) 
where N > 2 are also invariant under large gauge trans- 
formations. Such transformations are not continuously 
connected to the identity; rather, they generate infinitely 
degenerate perturbative vacua separated by action bar- 
riers that prevent classical transitions between them. 
Quantum mechanical tunneling can lead to a superposed 
ground state of these perturbative vacua, exponentially 
smaller in energy than any of them, since the tunnel- 
ing amplitude is itself suppressed exponentially by the 
height of the action barrier. A practical example is the 
QCD vacuum, where imaginary time (Euclidean space) 
solutions of minimum action (instantons) can be viewed 
as tunneling between adjacent vacua characterized by dif- 
ferent winding numbers n [14]. Taking into account this 
quantum tunneling, one usually writes the true vacuum 
state as a weighted superposition of all identical \n) vacua 

oo 

|6>>= Yl (!) 
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where for QCD, bounds from neutron electric dipole mo- 
ment studies [15] suggest that #qcd < 10~ 9 . However, 
for the purpose of illustration, let us assume that 8 is an 
unconstrained parameter of some SU(2)x gauge theory, 
which is broken at some high energy scale (we will gen- 
eralize our main results for arbitrary even Nf, where Nf 
is the number of fermionic doublets). Then, there are 
SU(2)x instantons, which represent tunneling solutions 
between the vacua with different winding numbers. 

Since this true ground state is lower in energy than 
any particular n-vacuum by only an exponentially small 
amount, the observed small but finite vacuum energy 
density can be explained if the universe has not yet set- 
tled down into the non-perturbative 0-state, but is still in 
one of the perturbative vacuum states. If one calculates 
the vacuum energy density in a 0-state or any n-state 
for pure gauge theory, t'Hooft's formula [16] for the one- 
instanton weight shows that the contribution of large size 
instantons diverges. The vacuum energy density in any 
\n) vacua can be expressed as [13] 

p v = (n\H\n) = 2Ke- Sa , (2) 

where So is the classical instanton action 8n 2 /g 2 and g is 
the coupling constant of the SU(2)x theory. Note that 
we assume the bare cosmological constant to vanish in 
the 6 — vacuum. 



In the dimensional regularization scheme [16], 

K = 2 . Vs -°/g ex P [--^ + f h ( W *) (3) 
+ 6.998435] . 

The integration over sizes R diverges but this can be 
tamed by introducing a physical cutoff scale po ~ M, set 
by an SU{2) scalar doublet <E> with a potential V($) = 
A(|$| 2 - M 2 /2) 2 , a la t'Hooft [16]. In this case, the pure 
gauge theory instanton should be replaced by the con- 
strained instanton solution [17], and the vacuum energy 
density in the n-vacuum is given by 

/87r\ 4 , 8- 2 
p v = (n\H\n) ~ M 4 e ^ . (4) 

If M is larger than Hi n f during inflation (for chaotic type 
inflation model H inf ~ 10 13 5 GeV [18]) , the inflation 
will naturally tend to wipe out inhomogeneities in $ over 
a Hubble volume. We will show that this argument con- 
strains the size of the instanton. 

Let us assume, for the purpose of illustration, 3 degen- 
erate vacua, |n = 0),|n = ±l). Note that we allow for the 
Universe to be in a state other than the non-perturbative 
|0)-state. The ground state of the Universe would then 
be 

|G) = -^(|0) + |1) + |-1», (5) 

and this would have a larger energy than if all winding 
number states \n = (—00,00)) were included (the en- 
ergy of the 9- vacuum is the least). Since the |0)-state is 
also an energy eigenstate unlike the |n)-state, we can, ex- 
pand each of the three \n) states in terms of the |0)-state. 
Given that (O\H\0) = p v (l - cos 6) *the probability dis- 
tribution function of the vacuum energy has a peak at 
2p v , corresponding to 6 = ir as shown in [13]. This 
peak can be interpreted as the value of the finite dark 
energy observed today. 

The inclusion of massive spin- 1/2 fermions in this 
SU(2)x gauge theory which admits instanton solutions, 
has interesting physical consequences. t'Hooft has shown 
that for fermion mass m <C 1/-R, where R is the typical 
instanton size, Nf species of spin 1/2 fermions contribute 
a factor [16] 

A f = (Rm) N i exp[~N f ln(MR) + 0.291746JV/] , (6) 



*This expression holds in the dilute gas approximation [19], 
corresponding to MR < 1, which is the regime we are working 
in (see below). 

*In QCD, 6 = 7r is a CP-conserving value, though here the 
is unrelated to strong interactions. 
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to the tunneling amplitude. Using the constrained in- 
stanton solution, the constant K now reads 



From the condition TH, 



K = A f 2 1 V g- s I 9f exp 



dR 



8tt 2 



-tt 2 R 2 M 2 (7) 
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+ 6.759189] . 



Let us suppose now that, as in [6], there are 2 fermionic 
doublets X and Y with different global quantum num- 
bers, that they acquire this mass through the sponta- 
neous breaking of the SU (2)x symmetry, and that the 
standard model quarks and leptons carry some quantum 
numbers under SU(2)x- Mixing between X and Y can 
occur at the quantum level via the SU(2)x anomaly, 
therefore, the quantum transition between the adjacent 
vacua must be accompanied by the decay process [6] 



X — > Y + SM quarks and leptons , 



(8) 



assuming that X fermions have a larger mass than Y 
fermions. The large mass of the X and Y (Kuzmin and 
Rubakov [6] assumed mx > 10 13 GeV) particles implies 
that the decay products (quarks and leptons) are highly 
energetic and can lead to the production of highly en- 
ergetic hadrons which can be the constituents of cosmic 
rays. 

The large mass value in our case implies the hierar- 
chy of scales R^ 1 > M > Hi n f. The first condition 
MR < 1 justifies the use of the M ^ instanton solu- 
tion, while the second condition M > H in f assures that 
gravitational effects on instantons are small during the 
inflationary phase and after inflation as well. Including 
the fermions, from Eq. (4), the vacuum energy in any 
perturbative |n)-vacuum is 



Pv 



m 2 M 2 e~ 



(9) 



where we have used the fact that MR ~ 1. As we have 
allowed for a small probability of transition T in the hori- 
zon volume, the condition THq 4 <~ 1 holds, where the 
tunneling rate per unit volume is 



^) m 2 M 2 e 

g 2 



(10) 



Using p v ~ 10" 120 M 4 and H 2 = p v /{3M 2 ) (where the 
subscript denotes the present era), we can solve for the 
coupling 



a 



4ir 



1 



44.1 



(11) 



Using this value of the coupling, we find from either of 
the equations Eq. (9) or Eq. (10) that, 
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1, we also have m 2 M 2 



a 4 Mp. We are proposing to explain the ultra- high en- 
ergy cosmic rays which can have energies upto the mass 
of the decaying AT-fermion, so it follows that one can 
expect cosmic rays of energies upto ~ 5 x 10 14 GeV 
(since M cannot be greater than M p ). We also note that 
M > Hi n f is satisfied. We may generalize our scenario 
to arbitrary number of fermion doublets Nf (which is re- 
quired to be even by the SU (2)x anomaly). Accordingly, 
Eq. (12) is modified to 



m N fM 4-N f _ 10 66 Gc y4 



(13) 



For example, for Nf — 4, m « 3 x 10 16 . For larger Nf, 
we find that the fermion mass scale is not considerably 
different. The fermion mass is indeed the prime result 
of this paper which clearly shows that the cascade decay 
of the fermions can give rise to ultra-high energy cosmic 
rays with energies greater than the GZK cut-off. 

Note that in the instanton mediated decay of fermions, 
the predicted mass turns out to be heavier than 10 13 GeV 
(for chaotic type inflation, the inflaton mass is around 
10 13 GeV in order to produce the right amplitude for 
the density perturbations and the spectrum [18]), and 
also greater than the observed spectrum from FLY's 
eye [21] and AG AS A [22]. Current analysis seems to 
be suggesting a relic fermion mass around 10 12 GeV 
[23]-10 14 6 GeV [24]. Particularly in our case, in order 
to excite the superheavy fermions at the very first in- 
stance, one has to rely on non-thermal production mech- 
anism for fermions after inflation [25] (see also the ap- 
pendix of Rcf. [26]). The upcoming experiment such as 
AUGER [27] will be able to see a considerable number of 
events above the GZK cut-off which will verify or falsify 
the energy scales which we predict here. If the inflaton 
coupling to the SU (2)x fermions is sufficiently large then 
the adequate abundance of such fermions will be the right 
candidate for cold dark matter. 

In summary, we argue that the problem of the cosmic 
dark energy, the small value of the cosmological constant, 
and ultra-high energy cosmic rays can have a common 
origin. We have shown that if the longevity of the X- 
particles is due to instanton-mediated decays, then the 
fermion mass, which sets the scale for the ultra-high en- 
ergy cosmic rays, must be larger than 5 x 10 14 GeV. Note 
that this is in accordance with the observed small cosmo- 
logical constant. 
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m 2 M 2 ~ 10 66 GeV 4 . 



(12) 
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